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Abstract 

The most important aim in tool path generation methods is to increase the 
machining efficiency by minimizing the total length of tool paths while the 
error is kept under a prescribed tolerance. This can be achieved by determin- 
ing the moving direction of the cutting tool such that the machined stripe is 
the widest. From a technical point of view it is recommended that the angle 
between the tool axis and the surface normal does not change too much along 
the tool path in order to ensure even abrasion of the tool. In this paper a 
mathematical method for tool path generation in 3-axis milling is presented, 
which considers these requirements by combining the features of isophotic 
curves and principal curvatures. It calculates the proposed moving direction 
of the tool at each point of the surface. The proposed direction depends on 
the measurement of the tool and on the curvature values of the surface. For 
triangulated surfaces a new discrete offset computation method is presented, 
which is suitable also for detecting tool collision with the target surface and 
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self intersection in the offset mesh. 
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1. Introduction 

This paper deals with three-axis milling, and presents computations on 
triangulated and analytical surfaces for determining the optimal moving di- 
rection of the tool at a given point of the surface. Basic requirements, such 
as the tool does not remove material from the target surface and the error re- 
mains under a prescribed tolerance, can be fulfilled by keeping the processed 
surface between the target (part) surface and its outer parallel (offset) surface 
at distance of a given tolerance. In the focus of this investigations are two 
problems of tool path planning; how to compute the offset of the mesh in a 
neighborhood of the actual point, and how to determine the optimal moving 
direction at any point of the surface represented by a triangular mesh. 

In a short survey mathematical approaches are mentioned dealing with 
geometrical factors in surface milling. The most frequently used tool path 
generation methods cut the surface by parallel driving planes in equal inter- 
vals or extract isoparametric curves in equally spaced parametric steps. If 
the surface has regions of different curvatures, both methods lead to uneven 
distances between tool paths, consequently, to fluctuating errors. The ma- 
chining error on the processed surface is measured by the height of the cusps 
(scallops) left between two adjacent tool paths (Fig. 1). Constant scallop 
height tool path generation methods for ball-end milling of analytical sur- 
faces are described in [29], [8] and [13], however, the required computations 
are expensive, and some technical questions (i.e. self-intersections) are not 
solved in all cases. For triangular meshes a constant scallop height method 
is presented in [TB] for computing the distances between slicing planes in 
three-axis milling with a ball-end cutter. 

Improvements of the frequently used isoparametric tool path generation 
method is suggested by including additional parameter curves in regions 
where the error is large [7], or by reparametrizing the surface in order to 
generate boundary-conform tool paths [12]. The boundary is followed by 
means of proximity maps in The application of two different tool path 
generation methods (using parallel driving planes or z-level contour parallel 
curves) is suggested by splitting the surface into "flat" and "steep" regions 
according to the angle between the tool axis and the surface normal in [BJ. 
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A method for approximating the processed part of the surface around a 
contact point, for generating the widest stripe along a tool path and for in- 
vestigating local and global millability with a given cutting tool is presented 
in [9] and (26]. The condition for collision- free manufacturing is expressed 
in terms of the curvature indicatrices of the part surface and the tool sur- 
face at their touching point. This local millability can be extended to global 
millability along an isophotic line. The angle between a given direction (the 
tool axis) and the surface normal does not change along an isophote, con- 
sequently the area of the contact surface of the tool is constant, and the 
abrasion of the tool is even. Isophotes are applied to partition the surface 
into regions, in which a generic side step (the distance between two neighbor- 
ing tool paths) can be computed from the prescribed maximal scallop height 
[5] . Such regions are computed for detection of interference between the sur- 
face and the tool (called usually gouging) in [28J. Isophote based tool path 
generation is proposed in [27] and an application to NC-machining in [10]. A 
numerical method for computation of isophotic curves on an analytical sur- 
face is presented in [TJ_]. Despite of the mentioned advantageous properties 
of isophotes these curves are not suitable milling paths due to their uneven 
density causing uncontrollable scallop heights. In general, each set of surface 
curves fulfills only some conditions for tool paths, but not all at the same 
time. 

Investigation of surface curvatures at the touching point of the target 
surface and the cutting tool concludes that the most efficient tool feed direc- 
tion in three-axis milling, which can minimize the total length of necessary 
tool paths, is the principal direction of the maximal normal curvature [2], 
[3] . Namely, this "steepest ascending" tool path leaves the widest track (pro- 
cessed stripe) on the machined surface, as also shown in [2S]. The same 
conclusion is formulated in [1] , where a torodial shaped cutter is investigated 
on a four and five axis machine. Vector fields are constructed for the optimal 
moving direction of the tool in [18J. A more flexible strategy for parallel plane 
milling considering the surface roughness is presented in |21j by introducing 
a fitness function proportional to the feed rate and transversal step. Then 
the acceptable machining directions are computed by optimizing the fitness 
function. 

In surface offsetting algorithms a difficult problem is to detect and to re- 
move self-intersections. For triangulated surfaces different offsetting methods 
have been published. Faces, edges and vertices are offset using half cylinders 
and spheres, then gaps are filled with blending methods and overlaps are 
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removed. By slicing this corrected offset triangular mesh with a series of 
planes three-axis tool paths are calculated [16]. By a local offsetting scheme 
facets are translated in the normal direction, convex edges to trimmed cylin- 
ders and convex vertices to trimmed spheres. Concave edges are removed 
from offsetting. Then in tool path generation the intersection with a driving 
plane results in circle arcs and line segments which are sorted and linked, 
overlapping portions are removed in order to generate a compound curve for 
cutting path [12] . [14"] . In a z- level contour tool path generation algorithm 
the triangulated surface is sectioned by a horizontal drive plane, then at each 
vertex of the section polygon an offset vector is created. A tool path is gen- 
erated sequentially along the section polygon while checking and removing 
intersections and using line fillets to connect offset segments pE]. Each tri- 
angle is offset in [20], the sliced segments in the drive plane are swept and 
chained into monotone polygonal tool paths. A different offsetting method 
is presented in this paper in Section 2. 

In Section 3 tool paths on triangulated surfaces are simulated by moving 
the tool in principal directions. In Section 4 proposals for optimal moving 
directions and their computation on analytical surfaces are presented. 

2. Mesh offsetting and computing the processed part of the surface 
at a point 

The result of the milling process with ball-end or toroidal tools is an 
approximation to the target surfaces, where the shape shows scallops (cusps). 
The maximal height of these scallops has to be kept within a given tolerance, 
i.e. the machined surface must lie between the target (part) surface and 
its offset by the given tolerance (Fig. 1). The piece of the tool surface 
lying in this region will be named contact surface of the tool. The curve of 
intersection of the tool surface and the offset surface is the boundary of this 
contact surface. The processed patch around a given touching point of the 
tool and the target surface is bounded by the projection of this boundary 
curve onto the part surface. The set of all processed patches should cover the 
entire part surface. In Fig. 2 a smooth surface, its offset, a touching ball-end 
tool and the processed patch on the surface are shown. 

In order to compute the processed patch at a given point, first the offset 
surface has to be generated. 

In the case, when the surface is represented by a triangular mesh, the 
offset will be computed in n normal planes. Each normal plane passes through 
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Figure 1: Milling with given tolerance 




Figure 2: The processed patch on the part surface 
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the touching point P of the tool, through the normal vector N of the actual 
triangle facet containing this point and an arbitrary direction in the plane of 
the triangle. The offsetting is carried out on the segments intersected by the 
normal plane from the triangle mesh (Fig. 3a). If the zth plane intersects 
the jth triangle in the segment Qij, Q 2 j, its offset is Qij + e'Nj, Q 2 j + e'Nj. 
Here e is the given tolerance, Nj is the normal vector of the jth triangle 
and e' = e cos a (a is the inclination angle between the normal vector Nj of 
the intersected jth facet and the ith section plane, e is shown by a thick 
segment and a is denoted by one arc in Fig. 3a). The position vector of 
a point is denoted by the corresponding boldface letter. In a convex region 
gaps arise between the offset segments, but in a concave region the adjacent 
segments intersect. Filling the gaps and removing self- intersections of the 
offset mesh are carried out by two-dimensional methods in each section plane. 
The obtained polygonal line intersects the normal section of the tool end in 
two points, Bu and B 2 j which are lying on the boundary curve of the contact 
surface of the tool and the offset mesh (Fig. 3b). The projections of B u and 
B 2 j onto the mesh give the endpoints of a diameter of the processed patch. 
This computation in the normal planes (i = 1, . . . ,n) results in 2n points 
lying on the boundary curve of the processed patch. In Fig. 4 the part of 
the offset surface within the ball-end of the tool is shown by 24 diameters 
above the base triangle. The boundary of the processed patch is drawn by a 
polygonal line on the mesh. Additional section planes can be included easily, 
if necessary. 

This local offsetting method is novel in the literature. It has the advantage 
that only two dimensional algorithms are necessary in the computation. Self 
intersection in the offset mesh and tool collision with the target surface can 
be detected by simple computation. Moreover, registering the directions of 
the normal planes, where such problems arise, a method for partitioning the 
surface could be developed for different milling tools and strategies. The 
applied edge-oriented polyhedral data structure defined on the mesh is very 
effective in determining plane sections of the mesh |22j . 

3. Moving the tool in principal direction 

First, the requirement is considered that the machined stripe along a tool 
path is the widest. Such tool paths will usually minimize the total length of 
milling paths [2], [3]. The processed patch by the tool around the touching 
point is determined by its n diameters lying on the surface. Now the tool will 
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(3a) Offset construction in a normal plane 



offset 




Figure 3: (3b) Computation of the intersection with the tool in a convex and a concave 
region 
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Figure 4: Part of the computed offset within the ball-end and the processed patch on the 
mesh of a quadratic surface. The actual triangle is lined. 

be moved in the direction perpendicular to the largest diameter. The next 
touching point of the tool will be chosen in this direction on the boundary 
of the processed patch. The suggested moving direction is similar in [26] in 
a different computational approach using Dupin-indicatrices. 

In Fig. 5 two tool paths on a quadratic saddle surface are shown with 
moving directions producing the widest stripe. The side step of the tool is 
chosen in such a way that the series of the processed patches overlap. 

Tool paths of the widest machined stripe are shown on a cylindrical sur- 
face in Fig. 6 and 7. A floating patch above one base triangle is also shown 
which is the part of the offset mesh within the ball-end of the tool. The mov- 
ing directions are shown by straight line segments. Each is emanating from 
the touching point on the actual triangle. The diameters of the processed 
patches are the shortest in the directions perpendicular to the generators in 
the case, when the tool touches the cylinder on the convex side, and they are 
parallel to the generators, when the tool is on the concave side of the surface. 
In the case of these "synthetic" meshes of a cylindrical surface the largest and 
shortest diameters of a processed patch are perpendicular to each other, and 
they are practically the principal directions. Investigations have been made 
on different synthetic meshes which have shown that the direction of the 
widest processed patch is very close to the principal direction of the largest 
curvature computed in a sufficiently large region around the contact point. 
For the estimation of the principal directions at a point of the triangulated 
surface a method is described in [23] and [23] by lying a circular disc onto the 
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Figure 5: Tool paths on the quadratic saddle surface shown in Fig 4 

mesh and computing the minimal and maximal chord lengths between the 
end points of its curved diagonals. In this concept the curvature values are 
ordered to the triangle facets of the mesh, and approximated surface normals 
are replaced by facet normals. 

The applied normal curvature estimation method constructed for trian- 
gular meshes is suitable for the computation on meshes also in such cases, 
where the known vertex oriented computations do not work. Namely, in the 
mesh of the shown cylinder all the vertices are situated on the surface bound- 
ary. To the contrary, frequently used vertex oriented methods require good 
vertices in the inside of the actual region. 

Even if the total tool path computed in the widest stripe direction can be 
the shortest, this kind of curves are not always appropriate tool paths due to 
technical reasons. Namely, the angle between the tool axis and the surface 
normal changes to much on a curved surface, which leads to uneven abrasion 
of the tool. The computed moving direction could be combined with other 
tool path generation algorithms as suggested also in [26J. 

4. Optimal moving direction and its computation on analytical sur- 
faces 

Now a technical requirement will be considered that the abrasion of the 
tool is even during the milling process. This requirement is fulfilled, if the 
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Figure 6: Tool paths on a cylindrical surface machined from the convex side 




Figure 7: Tool paths on the same cylindrical surface machined from the concave side 
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angle between the tool axis and the surface normal is constant along a tool 
path. Such curves are the isophotes on smooth surfaces. (The investigation 
of other material and kinematical factors is not the object of this paper.) 

As the isophotes are not suitable tool paths, a compromise will be sug- 
gested, and the tool will be moved from the contact point neither in the 
direction of the widest processed patch, nor along an isophote. A first idea 
was to choose the bisector of these two directions. The computation has 
been made on an analytical surface given by the function f(x,y), assuming 
that the piece of the surface to be machined is visible from the tool axis 
(i.e. the z-axis) direction [23]. In such cases the ^-projection method can be 
applied, when the processed patch is computed by projecting the boundary 
of the contact surface parallel to the ^-direction, and the step size is given in 
the xy-p\ane jlj. For the points B(x, y, f(x,y)) of the boundary curve of the 
contact surface the equation 



distance(C, B + eN(x, y)), N(x,y) 



(f' x (x,y)J' z (x,y),-l) 
yJfL 2 (x,y) + ff(x,y) + l 



yield, where r is the radius, C is the center point of the ball-end of the tool, 
and e is the offset distance. (, ) denotes the dot product. The points of the 
isophote passing through the point P satisfy the equation 

cos _1 (N p , a) = cos _1 (N(a;, y), a) , 

where a is the unit vector of the tool axis direction. In the specific coordinate 
system, where a\\z, this equation is equivalent to 

fL 2 ( x P>yp) + fy 2 (xp,yp) = f' x 2 {x,y) + fy 2 (x,y) . 

Both equations can be solved numerically resulting in a finite number of 
points lying on the required surface curves, i.e. on the boundary of the 
processed patch and on the isophote, respectively. 

Fig. 8 shows on a cubic surface how the processed patch around a contact 
point and the isophote through this point are situated. From the two common 
points of the curves two bisector directions can be determined (Fig. 9). One 
of them (which doesn't lead to a zig-zag curve) has been chosen for tool 
path direction. Few points of three generated tool paths with the processed 
patches are shown in Fig. 10. The step size is constant. The problem of 
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Figure 8: A processed patch around the contact point of the tool and the isophote through 
this point on a cubic surface 




Figure 9: Bisector directions hx and b 2 of the angle (W, P, Si) and (W, P, Sa), respectively. 
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Figure 10: Tool paths generated in the bisector of the widest stripe and the isophote 
directions. To each point the processed patch is shown. 

varying distances between the tool paths occurs here as well, as also in the 
frequently used methods ([5], [7]). 

A more reasonable idea for finding the optimal direction at any point 
of the surface is to consider the change of the angle between the tool axis 
and the surface normal (i.e. the inclination angle of the isophote passing 
through the actual point) while the tool is moving into a next position. This 
will determine the correction factor of the declination from the widest stripe 
direction. If the surface has small curvature, the surface normal doesn't 
change very much while moving in the principal direction of the maximal 
curvature. Therefore, the tool can be kept on this surface curve, because the 
abrasion will be even. In the case of larger curvature values the tool should 
be moved rather towards the corresponding isophote direction. 

In a generic case the suggested moving direction will be determined be- 
tween the widest stripe direction and the direction pointing to the isophote 
point on the boundary of the processed patch c in such a way that it divides 
the arc between these two directions in the ratio (1 — cos/3) : cos/3. /3 denotes 
the difference of the inclination angles of the surface normals to the tool axis 
at the two corresponding points on the processed patch boundary. In Fig. 
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11 the vector w is pointing into the direction of the widest stripe and s is 
pointing to the point S of the isophote passing through the contact point P. 
The boundary c is represented by 2n points Sj which are the projections of 
the points lying on the boundary of the contact surface denoted by and 
B 2 j in Fig. 3b (n is the number of the section planes). The isophote point S 
on c can be approximated by one of the points Sj, in which for the angles 

(N Si ,a)Z«(N P ,a)Z 

yields, computing from the condition 

| cos _1 (N s , a) — cos _1 (Np, a)| = min\ cos _1 (N Si , a) — cos -1 (Np, a)| , 

Sj e c, i = 1, . . . , 2n , 

where Ng. is the normal vector of the triangle containing the point That 
means, the direction of Ng should fit on the cone determined by the isophote's 
angle of Np (denoted by single arc in Fig. 11). The vectors Ng, Np, Ng. 
and a are unit vectors. 

Let (3 denote the difference of the angles which are formed by the surface 
normals and Ng to the axis direction a. 

P = |cos -1 (N w ,a) -cos -1 (N s ,a)|, ||N W || = 1, ||N S || = 1, ||a|| = 1. 

The proposed new moving direction q is determined by the arc length of the 
processed patch boundary between w and s: 

arc(q, w) = (1 — cos 0) ■ arc(s, w). 

Consequently, a zero, or a very small (3 does not change the moving direction 
w. This P depends also on the radius of the ball-end. The bigger the tool, the 
bigger is the processed patch. Therefore, the change of the inclination angle 
of the normals is larger. We remark that the factors included in the fitness 
function introduced in [21 J have the same effect on the moving strategy. 

This path correction does not influence the abrasion of a ball-end tool, 
it has importance in milling with other cutters. The computation for a 
toroidal cutter does not differ from the presented computation above, but 
the demonstrating figures have been made with ball-end tools. 

Fig. 12 shows a tool path computed by the suggested method. The path 
in the widest stripe direction (practically a parallel circle of the cylinder) 
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Figure 11: The proposed moving direction q is computed between the widest stripe direc- 
tion w and the isophotic direction s. 

is shown with the processed patches around the computed contact points. 
The corrected tool path situated lower shows a declination. The radius of 
the ball-end of the tool is in this case eight times the average size of the 
triangles. The declination of the new path is less with a smaller radius. In 
the case, when the radius is only four times the average size of the triangles, 
the tool will move approximately along the path of the widest stripe. The 
tolerance was one third of the tool radius in this case, which is the distance 
of the floating piece of the offset and the part surface (see on the left-hand 
side in the figure). Note that the tool paths computed on a cylinder surface 
in a completely different way in [21] show a similar shape. 

The next examples are computed on parametric surfaces. Though the 
surface is presented by a vector function r(u,v), (u,v) G [wi, 1 ^] x [^1,^2], 
which is at least twice differentiable in a sufficiently large neighborhood of 
P, the computation is carried out in n normal planes producing 2n points 
of the surface curve c bounding the processed patch. In this way, the time 
consuming and unstable process of computing the projection of a point onto 
the surface will be replaced by a very effective method. 

Let's say, the actual tangent direction vector is Sj, i G {1, ...n}. The 
normal plain containing the required point Sj is determined by the vectors Sj 
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Figure 12: Tool paths on the cylinder in the widest stripe direction and in the proposed 
moving direction. 

and the surface normal Np. The decomposition of Sj in the tangent plane is 



Let the point to be projected onto the surface be Bj lying in this normal 
plane. First, and 6j are computed as solutions of the system of linear 
equations from the decomposition of Sj. Then series of points are generated 
in the parameter domain with an appropriate At, which are 



The corresponding points on the surface are determined by the vectors r(uk, Vk) 
The required point Sj is determined by the parameter values (uk, v^) for which 
the distance of the surface point r(uk, i>k) to Bj is minimal (k e {1, . . . , K}. 

Numerical tests have shown that the results of the presented approxima- 
tion of the projection onto the surface are satisfying compared to the "real" 
projections computed by the program package Mathematica by searching for 
a nearest surface point to the one to be projected without any constraints. 
The difference between the results of the two methods is within the error 
bounds of the computation in other procedures of the complete task. 
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Figure 13: Tool path on the torus shown with the processed patches around the computed 
points. 

In the Figures 13 and 14 the results of this computation on a torus are 
shown. The smaller patches are generated with a smaller tool radius. Here 
the tool path deviates little from the parameter curve of the meridian circle, 
which is the path of widest stripe with tangents in the principal directions 
of the greater principal curvatures. A larger tool radius shows a significantly 
bigger deviation. 

The results in the examples correspond to expected ones, the direction 
choice optimizes the tool inclination angle for maximal material removal rate 
in each cutter contact point. The suggested locally optimal moving direction 
is depending as well on the surface curvatures at the actual point, as on the 
tool's size. 

In the following analysis three examples are compared. The table shows 
the intervals of the Gaussian curvatures in the processed region, the change 
of the tool inclination angle in the steepest ascending (i.e. maximal curva- 
ture) direction and in the proposed direction, respectively, in one step. The 
last column shows the decreasing rate of the width of the processed patch 



17 



Figure 14: Tool path on the same torus with larger tool radius. 
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compared to that of the widest stripe. The three surfaces are a torus (Figures 
13 and 14), a saddle surface and a trigonometric surface. Figure 15a shows 
the path of a relatively large tool compared to the principal curvatures on 
a saddle surface. Therefore, it deviates very much from the deepest ascend- 
ing direction (Figure 15b), and goes closed to the isophotic curve (Figure 
15c). The bumpy trigonometric surface in Figures 16a and 16b shows that 
the computed directions are stable also in more curved regions. Here the 
steepest ascending direction crosses the bump (not shown), and an isophotic 
curve is shown in Figure 16c. The proposed tool path goes between them, 
because the tool's size compared to the principal curvatures is smaller then 
in the case of the saddle surface. 
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5. Conclusions 

Geometric considerations about generating milling tool paths have been 
presented. Local moving direction of the tool has been computed, when 
the processed stripe is the widest, and the change of the inclination angle 
of the tool axis is minimal. A correction factor to the steepest ascending 
direction has been determined from the change of the isophotic angle of the 
surface normal with respect to the tool axis while moving the tool from a 
given point into the next position. This strategy provides to develop milling 
strategy with optimal tool paths, when only these two geometric requirements 
are considered. The presented new local offsetting method for triangular 
meshes solves basic problems, as computing the processed patch boundary 
and detecting self intersection or gauging. For the solution of different other 
technical problems further investigations are necessary. 

On analytical surfaces the computations and the figures have been made 
by the symbolic algebraic program package Mathematica. The program de- 
veloped by the author for the presented computations on discrete surfaces is 
not for professional applications. A lot of numerical problems may arise in 
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Figure 16: Tool paths in a more and less curved region and an isophote on a trigonometric 
surface. 



the computation with real triangular meshes. The solution of such problems 
was not the task of this work. 
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